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Abstract 

In this paper, we perform a complete one-loop computation of the short-distance Wilson 
coefficients for B s c i — B s ^ mixings and B s d £ + £~ decays in the Manohar-Wise model, 
which extends the SM scalar sector by a colour-octet and weak-doublet scalar. Based 
on these calculations, combined constraints on the model parameters are derived from 
the current flavour data, including the B S( i — B sc i mixings, B S( i —>• B — > X s j, 

B —> K * 7 , B —> fry, and Z —> bb decays. The future sensitivity to the model is also 
explored in the observables achievable with 50 fb _1 of LHCb and 50 ab ^ 1 of Belle-II 
data. We find that the Manohar-Wise model could explain the current data, especially 
when the couplings of the charged colour-octet scalars to quarks are complex, with the 
resulting Xmin being significantly smaller than that of the SM. Finally, we investigate cor¬ 
relations among the two isospin asymmetries A(iL* 7 ) and A(/ry), as well as the averaged 
time-integrated branching ratios B(B S ^ some of them are found to be quite 

strong and could provide, therefore, further insights into the model, once more precise 
experimental measurements and theoretical predictions for these observables are available 
in the future. 
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1 Introduction 


The discovery of a new boson, with a mass close to 125 GeV, by the ATLAS (T] and CMS |2 
collaborations stands as a remarkable success of the Standard Model (SM) of particle physics. 
The measured properties of this boson are so far in agreement with those of the SM Higgs |3 ^5], 
suggesting that the electroweak symmetry breaking (EWSB) is probably realized in the most 
elegant and simple way, i.e., via the Higgs mechanism implemented through one scalar doublet. 
It is, however, noted that none of the fundamental principles of the SM forbids the possibility 
of an enlarged scalar sector associated with the EWSB. A natural question we are now facing 
is then whether the discovered state corresponds to the unique Higgs boson predicted by the 
SM, or it is just the first signal of a much richer scenario of EWSB. 

Among the many possible scenarios for new physics (NP) beyond the SM, the two-Higgs- 
doublet model (2HDM) |6] provides a minimal extension of the SM, by adding a second scalar 
doublet to the SM field content. It can easily accommodate the electroweak (EW) precision 
tests and give rise, at the same time, to a very rich phenomenology |7|. As a potential theory 
of nature, the 2HDM is very interesting on its own, since it allows for CP violation beyond 
what is provided by the SM. It is also helpful to gain further insights into the scalar sector of 
supersymmetry and other models that contain similar scalar contents. The direct search for 
the scalar spectrum of the model at high-energy collisions or through indirect constraints via 
precision flavour experiments are, therefore, an important task for the next years. 

Within the SM, the flavour-changing neutral currents (FCNCs) are forbidden at tree level 
and, due to the Glashow-Iliopoulos-Maiani (GIM) mechanism |8|, are highly suppressed at loop 
level. In the most general version of 2HDM, however, unwanted FCNCs appear even at tree 
level, which represents a major shortcoming of the model. The hypothesis of natural flavour 
conservation (NFC) is the usual way out to this issue. By limiting the number of scalar doublets 
coupling to a given type of right-handed fermion to be at most one, the absence of dangerous 
FCNCs is guaranteed ||9|. This can be explicitly implemented via a discrete Z 2 symmetry acting 
differently on the two scalar doublets, leading to four types of 2HDM (usually named as type-I, 
II, X and Y models), which has been studied extensively for many years |7|. 

Another efficient way to guarantee the smallness of FCNCs is to impose the principle of 


minimal flavour violation (MFV) (10 12 , a concept that can be traced back to Refs. 13,14 and 
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has become very popular during the past decade. It amounts to assuming that all the flavour- 
violating interactions, including those mediated by the electrically neutral scalars, are controlled 


by the Cabibbo-Kobayashi-Maskawa (CKM) matrix 15 16 , as happens in the SM. This can 


be implemented by requiring all the scalar Yukawa couplings be composed of the SM ones 


Y u and Y D . It has been shown in Refs. 17,18 that, under the MFV hypothesis, the allowed 


SU(3)c®SU ( 2 )(l)^- representations of the second scalar that couples to quarks via Yukawa 
interactions are fixed to be either (1, 2)i/ 2 or (8, 2)i/ 2 ] namely, the second scalar can be either 
colour-singlet or colour-octet. Examples of the former include the aligned 2HDM (A2HDM) [T9] 
and the four types of 2HDM reviewed in Ref. |7|. In the following, we shall refer to the 


2HDM with the second scalar being colour-octet as the Manohar-Wise (MW) model 17 . A 


characteristic feature of the MW model is that its scalar spectrum contains, besides a CP-even 
and colour-singlet Higgs boson (the usual SM one), four colour-octet particles (one CP-even, one 


CP-odd and two electrically charged), giving rise to many interesting phenomena 17,20-25 


I 11 this paper, motivated by the latest experimental data on B s ^ — B Std mixings and B S) d —> 


fi + fi decays 26, 27 , we shall perform a detailed study of these processes within the MW 


model. While the B sd — B s d mixings have already been addressed within the model 17 , it 
should be noted that the Wilson coefficient (Eq. (26) in [l7|) is obtained only in the limit of 
heavy charged colour-octet scalars and with zero external momenta. In a more general case 
for the model parameters, contributions from the other operators become non-negligible and, 
in order to get a gauge-independent result, the box diagrams should be calculated by keeping 
the external momenta up to the second order. The same observation is also applied to the 


rare leptonic B sd —> £ + £ decays, as demonstrated in Ref. 28 . The current experimental data 
on the averaged time-integrated branching ratios of B s d —> /i + //“ decays, averaged over the 


CMS 29 and LHCb |30| measurements, read 31 


6(5 s ^/iV) = (2.8+°;^) x 10 


-1-9 


B(B d -> /iV) = (3-9±i;l) x 10 


,-10 


( 1 . 1 ) 


which are in remarkable agreement with the latest updated predictions within the SM 32 


B(B S fj+fj,-) = (3.65 ± 0.23) x 10 ~ 9 , B{B d -)• ^ n~) = (1.06 ± 0.09) x 10" 1U . (1.2) 


-10 


All the experimental and theoretical progresses will lead to more stringent constraints on physics 
beyond the SM. Combining these processes with the interesting B —>• W. 7 , B —s- K* 7 , B —>■ 
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P7) and Z —>• bb decays 21 22 , we shall investigate the combined constraints on the model 


parameters from the current flavour data. In addition, the future sensitivity to the model is 


also explored in the observables achievable with 50 fb 1 of LHCb 


33 


and 50 ab -1 of Belle-II 

data [34] (the so-called “stage II” projection introduced in Ref. [35]). 

The outline of this paper is as follows. In Sec. [2] we give a brief review of the MW model. In 
Sec. [3] we summarize the theoretical framework used to calculate the flavour observables both 
within the SM and in the MW model. A complete one-loop computation of the relevant Wilson 
coefficients for B s d — B sd mixings and B sd —* l + i~ decays in the MW model is also described 
in this section. Detailed numerical results and discussion are then presented in Sec. [4] Our 
conclusions are finally made in Sec. [5] Explicit analytical results for the Wilson coefficients 
relevant to B sd — B S)d mixings are given in Appendix |Aj 


2 The Manohar-Wise model 


In the MW model 17 , the scalar sector of the SM is supplemented with a colour-octet and 


weak-doublet scalar. The two scalar fields in the model can be parametrized by 


H = 


co 


S A = 


( 2 . 1 ) 


±(v + h ~i z °)J VA ( S £ + iSf) y 

with v = (y/^Gp)- 1 / 2 ~ 246 GeV. Here H is identified as the SM scalar doublet, with uA and 
z° being the three would-be Goldstone bosons, whereas h is the SM Higgs. S A denotes the 
colour-octet and weak-doublet scalar, with A being an adjoint colour index; S±, S Rj and S A 
denote the electrically charged, neutral CP-even, and CP-odd colour-octet scalars, respectively. 


2.1 Yukawa couplings to fermions 


The Yukawa interactions of the two scalar fields with the SM fermions are given by 17,22 


Cy = Ql(Y d H + Y d S A T A )d° R + Ql(Y u H + Y u S A T A )u° R + L° L Y l He R + h.c., 


SO (\ru 


'ArpA\ n ,0 


z 0xr£ } 


( 2 . 2 ) 


where H = and S A = ia 2 S A *, with 02 being the Pauli matrix. Q R and L° l denote the left- 

handed quark and lepton doublets, and u Rj d R and e R the right-handed up-type quark, down- 
type quark and lepton singlets, respectively, in the weak interaction basis. The fundamental 
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representation in colour space of the SU(3)c generators T A {A = 1, • • • , 8) acts on the quark 
fields, and determines the colour nature of the second scalar doublet. 


All fermionic fields in Eq. (2.2) are written as 3-dimensional flavour vectors and the Yukawa 


couplings Y* and Y? (/ = u, d, £) are, therefore, general 3x3 complex matrices in flavour space. 
According to the MFV hypothesis, the matrices Y u,d should be composed of the combination of 
Y u and Y d , and transform under the SU(3)q l <g) SU(3)u R <8> SU(3)d r flavour symmetry in the 


same way as Y u,d themselves 17 . This can be achieved by requiring the alignment in flavour 


space of the Yukawa matrices 19 


Y d = r) d Y d , Y u = rj* u Y u , (2.3) 

where the two parameters r]f (/ = u, d ) are generally arbitrary complex numbers. The align¬ 


ment condition also guarantees automatically the absence of tree-level FCNCs 19,36,37|. Ap 


plying the SM unitary transformations to rotate the fermionic fields from the interaction to the 
mass-eigenstate basis, one can finally obtain the Yukawa interaction of physical colour-octet 
scalars with quarks in the mass-eigenstate basis [l7| 


C Y = - ^ u u l ^T A u l L S A - V2rj d (P L T A — d i R S ( 


m n 


v 


v 


^u^T^diSi - V2r ]d u i L 1 ^T A Vi j d : ’ R S A + h.c. 


m a 

V 


(2.4) 


where S A = Sr and V.\j (i — u,c,t and j = d, s, b ) is the CKM matrix element. 

2.2 Scalar self-couplings and couplings to gauge bosons 

The most general and renormalizable scalar potential of the MW model is given by |17| 

\ / 2 \ ^ 

V =- ( H^Hi - V — ) + 2m|TrA t h5' l + S 3 + X 2 H^H j TrS^ j S i 

+ hsH^H^TvSiSj + X^TrS^SjSi + A,,// r lV.S' T '.S',,S; ; + h.c. 

+ A 6 Tr S^SiS^Sj + XrTrS^SjS^Si + AgTr^^Tr^^- + XgTvS^SjTrS^ S t 

+ XwTrSiSjTr&'tf* + A n Tr (2.5) 

where the SU{2) L indices on the scalar doublets have been explicitly displayed. Traces are 
taken over the colour indices and the notation S = S A T A is used. Due to the Hermiticity of 
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the scalar potential, all the parameters are real except for A 4 and A 5 (A 3 has been made real by 
a phase rotation of the S fields) [l7|. From this potential, one can easily obtain the trilinear 
coupling involving two charged colour-octet scalars and the SM Higgs boson h 

9sfs B h — ~ 2 ^ AB - (2-6) 

which is the only cubic coupling relevant in our calculation. 

The scalar couplings to the gauge bosons arise from the gauge-kinetic Lagrangian 

C kin = {D^H) ] D»H + {D,S A y D^S a , (2.7) 


where D M denotes the SM gauge-covariant derivative. From this Lagrangian, one can easily get 
the trilinear coupling of two charged colour-octet scalars with the Z boson 

cos 2 9 W 


9sfs B zv —9 0 n 
+ - 2 cos 0 W 


P^B PgA ) $ABi 


( 2 . 8 ) 


where g is the SU(2) L gauge coupling constant, and 6w the weak mixing angle. p^ A and p^ B 
denote, respectively, the incoming momenta of the two charged colour-octet scalars. 

As all the other couplings are irrelevant to our calculation, they are not mentioned here. It 
is, however, noted that all the Feynman rules in the MW model have been cross-checked with 


the package FeynRules 38 


3 Theoretical framework for flavour observables 

3.1 B s , d - B s , d mixings 

The B-meson mixings are governed by the SM box diagrams with exchanges of up-type quarks 
and W ± bosons, and the box diagrams mediated by up-type quarks and charged colour-octet 
scalars in the MW model, as shown in Fig. [TJ The resulting effective weak Hamiltonian reads 

«la B|=2 = ^ m 2 „, V CMQ, + h.c„ (3,1) 

i 

Here Gf is the Fermi coupling constant, and Ci(p) = l / c KM C'j(/i) with Vq KM being the CKM 
factor and Ci(p) the scale-dependent Wilson coefficients. When QCD renormalization group 
effects are taken into account, there are totally 8 operators responsible for the B s ^ — B sc i 


39 
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Figure 1: Box diagrams for the B q — B q mixing in the unitary gauge both within the SM (the first 
one)and in the MW model (the last three). Crossed diagrams have also been taken into account. 


mixings, which can be split further into 5 separate sectors according to the chirality of the 


quark fields they contain (q = d or s ) 139 40 


Q\ ll = (g°7 ' l >PLb a )(<?'fP L bP), 

Qi R = (q a ^P L b a )(frPRb P ), Q RR = {?P L b a ){tfP R bP), 
Q\ LL = (fW)W), Ql LL = (q^PLb^a^P^), 


(3.2) 


where a,/3 are the colour indices, a /w = ^[ 7 ^, 7 ^] and Pl,r — |(1 T 7 . 5 )- The remaining two 
sectors (Q\ RR and Qf RR ) are obtained from Q\ LL and Qf LL by interchanging P F and Pr. 
Within the SM, only the operator Q\ LL contributes and the corresponding Wilson coefficient 


up to next-to-leading order (NLO) is given by 41 42 


CyLLil^w) = 4A 2 t Sww{xt) + ~ , i '' ~ D S M{x t ,x 


VW J 


(3.3) 


with A* = V t bV t * q , x t = and x = jj^/rri^. The leading order (LO) coefficient Sww(%t) 


is the known Inami-Lim function 43 , and the NLO coefficient DsAt(xt, x^w) reads 42 
Ds M (x t , Xfi W ) = C F \ L (1 ' SM \x t ) + 


61n(aw) ( xt-^- ) +3 


Sww(xt) 


+ C A \ L^ SM \x t ) + 


6 In (xpw) + 5 


S W w{x t ) c , 


(3.4) 


where C F = 4/3 and Ca = 1/3. Explicit expressions for S ww (x t ), L^' SM \x t ) and L^' SM \x t ) 


can be found in Refs. 41,42 . Note that the SM box diagrams are evaluated with zero external 


momenta, which is a good approximation for B s ^ — B S) d mixings within the SM 


44 


The additional contribution to the effective weak Hamiltonian arises from the SM box 
diagrams with the W ± bosons replaced by the charged colour-octet scalars (the last three ones 
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in Fig. [Tj). These diagrams depend on the two parameters rj u and r/^, which characterise the 
Yukawa couplings of charged colour-octet scalars to quarks. For most general values of these 


parameters, especially when r/d/r] u ~ rrit/mb , each term of the second line in Eq. (2.4) can give 


a comparable contribution and should be, therefore, taken into account simultaneously. In such 
a specific case, in order to obtain a gauge-independent result, it is necessary to keep at least 


the external 6-quark momenta up to the second order 28 . This is contrary to the SM case 


where all the momenta of external quarks can be safely set to zero. 
Following the same computational procedure as used in Ref. 


, we set the light-quark 


masses md )S to zero; while for the 6-quark mass mj, we keep it up to the second order. As 
the external momenta are much smaller than the masses of internal propagators, the Feynman 


integrands are firstly expanded in external momenta before performing the loop integration 45 

l 2 + 2(k-l) 


4 (k-iy 

(k + l) 2 — M 2 k 2 — M 2 f k 2 — M 2 ' (k 2 - M 2 ) 2 


1 - 


+ 


+ 0(1*/M 4 


(3.5) 


where M denotes a heavy mass, k is the loop momentum and l an arbitrary external momentum. 


After factorizing out the external momenta in Eq. (3.5), we are left with integrals that depend 


only on the loop momentum k and the heavy masses m,> They can be further reduced, with 


the help of partial fraction decomposition 46 

1 


(q 2 — m\)(q 2 — m 2 ) m\ — m\ L 


(3.6) 


q* — q * — rn, 2 

to the ones in which only a single mass occurs in the propagator denominators. Finally, after 


the tensor reduction 47 , the only non-vanishing one-loop scalar integrals take the form 48 

n—D/2 


d D k 


(— 1 ) n i r(n — D/2) ( 1 


77T 


(3.7) 


J ( 2n) D ( k 2 — m 2 ) n (4(/r) D / 2 T(n) 

with an arbitrary integer power n and with m ^ 0. In addition, the naive dimensional regular¬ 
ization scheme is employed to regularize the divergences appearing in Feynman integrals. 

The LO non-vanishing Wilson coefficients in the MW model at the matching scale (j,$ ~ 
0(7ti s a), where heavy degrees of freedom are integrated out, are given by 


CdEM = A JcSV) + 2A C A + ttW, 

<Y5ufe) = A?CS>s) + 2\ c X t C^gM + A 


(3.8) 


(3.9) 


























(3.10) 


PS NP 
^SRR, 2 


i^s) 


\2p/NP,tt 
A t ^SRR ,2 


(/^s) + 2X c X t C 


NP,ct 
SRR, 2 


M + A^> s ), 


where A c = V cb V* q and the explicit expressions for Cy^, Cg^\ and are listed in 

Appendix [A} Here we have used the Fierz identities 


(q^P L T A b)(qrPLT A b) = i Q\ LL , 

0 qi^P L T A T B b)(qYP L T B T A b ) = iigP 1 , 
(qP B T A b)(qP„T A b) - -y Qf«K + I Qf™ 


(3.11) 

(3.12) 

(3.13) 


It should be noted that the limit rn u —> 0 and the unitarity relation X u + A c + X t = 0 have 
been implicitly exploited during the calculation. To make sure the gauge independence of our 
results, we have performed the calculation both in the Feynman and in the unitary gauge. 
Together with the analytical formulae for the QCD renormalization group factors given in 


Ref. 39 , one can obtain the corresponding Wilson coefficients at the lower scale ^b ~ 0(m b ). 
The effective weak Hamiltonian for I AH I = 2 transitions in the MW model reads 


n 


|AB|=2 

eff 


G 2 F 

167r 2 


771 


w 


c?lM + Q\ LL + cg&MQt 


VLL 


}SRR 


+ GsRR,2 {P>b)Q2 RR | + h- c - ; 


(3.14) 


from which the off-diagonal matrix element for B q — B q mixing is derived as 

c™M + c»ZM] (B q \QX LL ^b)\B q ) 


(B q \H^ Bl ~ 2 \B q ) = 


167T 2 


W 


+ CZrM) <B,IQ?*'V)|B,> + (B,|Qf'"‘( W )|B 5 > \ , (3.15) 


SRR/ 


with the operator matrix elements (B q \Qi(/j, b )\B q ) given, respectively, by 39 


(B,\Q\ LL (q t )\B q ) = 1 ''''(W). 


(B«IQf RH (w)|B,> = -^7 


m Ba 


24 \m b (fi b ) + m ? (// 6 

(b«\QI rr M\b«) = -1 ( 


m B q f Bq 5f K "(/i b ), 


SRR/ 




m b (n b ) + m q (ii b ) 


m Bq f B B 2 RR {iXb)- 


J q J ttq 


(3.16) 

(3.17) 

(3.18) 
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Here fp is the Z^-rneson decay constant, and B\ LL , Bf RR and are the non-pertnrbative 


Bag parameters, whose values can be obtained in lattice QCD 49 . The mass difference AM B o 
is then given by 


AMgo = 2 




(3.19) 


which can be used to constrain the model parameters, taking into account the current precise 


experimental measurements 26,27 . 


3.2 B s ,d —>■ decays 

The most general effective weak Hamiltonian for the purely leptonic B s d —> l + l~ decays can 


be written as 50 


G F m W 


7T 


10, S,P 

v «y» E (c i o i +c;o')+h. 


c. 


(3.20) 


with the corresponding operators given by 

Oiq = (Ql^Lb) (h p ld) , 
m e m b 


O s = 


M 2 W 


(qPnb) {U ), 


° p = IP ^ ’ 

lvl w 


O'lc = (<I1,PrV) , 

mi m e m b . _ . 

°S = -JJ2~ \d P L b ) (tz) , 

m w 

o' P = ^ mb) (hd) , 

m w 


(3.21) 


where t = e, /i, t and q — d,s. In this paper, we shall neglect contributions from the operators 
0[, because they only give contributions proportional to the light-quark mass m q . These rare 
processes are then governed only by the operators O i 0 , Os and Op in our approximation. 

Within the SM, the dominant contributions to the rare B s ^ —» i + d decays originate from 
the W-box and Z-penguin diagrams shown in Fig. [2j which generate the Wilson coefficient 


where 


Y 0 (art) = ? 


n ? CD Y 0 (x t ), 

(3.22) 

x t 4 3x t 

+ In x t 

x t -1 [x t -1 y 

(3.23) 


is the one-loop function calculated firstly in Ref. 1431. The factor r^y W accounts for both the NLO 


EW matching corrections 51 and the logarithmically enhanced QED corrections originating 
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Figure 2: SM Feynman diagrams contributing to the rare B q —> /x + /r decays in the unitary gauge. 


from the renormalization group evolution 32,52 , while ?/y CD stands for the NLO 


53 


54 


and 


next-to-next-to-lcading order (NNLO) [52 QCD corrections. 


As detailed in Ref. 28 , the W -box and Z-penguin diagrams also give rise to the pseudoscalar 


coefficient Cp. On the other hand, the scalar coefficient Cg is generated exclusively from the 
W -box and Higgs-penguin diagrams shown in Fig. [2j The complete SM expressions for Cp and 
C's in the unitary gauge are given, respectively, by |28] 


/iSM _ /obox, SM | 

^ P ^P, unit ary ' ^P, unitary 


yZ penguin, SM 


l 

24 

.2 

’W 


ay (36a; 3 — 203a; 2 + 352ay 


6 (ay - l) 3 


209) 17xj - 34a: 3 + 4a: 2 + 23ay 


{x t - l) 4 


ay(18a; 3 - 139a: 2 + 274ay - 129) 24a; 4 - 33a: 3 - 45a; 2 + 50ay - 8 


36 

ySM _ ^ybox, SM 

mitar 

3 x t 


/ObM ^yDOX, aivi | ^-y 

S ^ S, unitary ' ^ S, unitary 


2{x t - l) 3 + 

h penguin, SM 
S , unitary 

x t (x t + 1) (ay — 2) (3a: 2 — 3ay + 1) 


(x t ~ l ) 4 


In x t 


In ay 


(3.24) 


In x t , 


(3.25) 


8 a; hsM 48(ay - l) 2 24(ay - l) 3 

where xj lSM = and m^o is the mass of the SM Higgs boson. 

In the MW model, the colour-octet scalars do not couple to the leptons and there is, there¬ 
fore, neither additional box diagrams with charged colour-octet scalar exchanges nor penguin 
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diagrams with neutral colour-octet scalar exchanges. The only new contributions to the rare 
B S) d — > £ + £~ decays come from the penguin diagrams with the Z boson and the SM Higgs boson 
exchanges, as shown in Fig. [3] The new Z-penguin diagrams involve the charged colour-octet 
scalar exchanges, and provide additional contributions to the Wilson coefficients C\q and Cp\ 


/-yZ penguin, NP 
^10, unit ary 


= \v u n 


+ 


Xs 


s-yZ penguin, NP _ 

^P, unit ary 


6 [x s - x t ( Xs - x t ) 
Xt 


(In xt — In xs) 


(3.26) 


3 (x s -x t y 
2 


VdVu 


Xt + X S , x t x s n . , 

-v-h 7-r (lnxs ~ mx t ) 

2 ( x s ~ x t ) 


I Vu 


6 (x s -x t ) 

Zs^xt 
9(^ s - Xt y 


4 - 8 x s x t - I7x 2 t x 2 t (3x5 + x t ) , , ' 

H - 7 - 7 — (hi Xs - m x t ) 


VdVn 


6 (xs - Xt) 

5 x t - 3x s , x s (2 x s - 3x t ) 


+ 


(xs ~ x t ) 


(lnxs — lnx t ) 


+ 


\Vv 


I 2 

1 

17x| — 64 xsx t + 71x 2 

4x| - 12 x| x t + 9x s x 2 + 3xl nn r ln T i ) 

- Xt) 

6 

msj m x t ) > 
(x s - x t ) J 


+ \Vu\ 2 (1 - s 2 w) 


Xt 


x s (lnx '5 - lnx t ) + x t - x s 


(3.27) 


3(x s - x t ) _ 

where xs = m 2 sA /m ^ and rn S A is the mass of the charged colour-octet scalars. While the 
new SM Higgs-penguin diagrams involve also the charged colour-octet scalar exchanges, they 
contribute only to the scalar Wilson coefficient Cs- 


penguin, NP v 2 X t \\ 
^unitary “ m 2 g (^ _ Xf ) 


+ 1 ^ 


+ 


V 

12 

2 x+ 


Wu 


x t 


_(x s -x t ) 


Xt 


(In xs — lnx t ) — 1 


xs ~ 3x t 


-7T (I 1 IX 5 — 111 X t ) + „ , 

[2(x s -Xt) 2K ’ 4 (x s -x t ) 


3x 

+1 n 


h-SM 

2 


VdVu 


Xt 


Xt 


Xs 


(x S - Xt) [ (x S - Xt) 


2 (x s -x t y 


Xs + X t 


x s x t 


(x S - x t ) 


(I11X5 — lnx t ) 


(In xs — In xt) 


(3.28) 


With Eqs. (3.22)—(3.28), the final Wilson coefficients in the MW model can be written as 

C 10 ~- 


(SM , ^-yZ penguin, NP 

10 "r '-'10, unitary J 

(3.29) 

(SM , sih penguin, NP 

S ' ^5, unitary J 

(3.30) 

(SM 1 rPL penguin, NP 

P ' P, unitary 

(3.31) 
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Figure 3: Penguin diagrams contributing to the rare B q —> /P/x decay within the MW model. 


With the operators defined by Eq. (3.21), none of the above three Wilson coefficients is affected 
by the QCD renormalization group evolution. Together with the hadronic matrix elements, 


(%7/W5&|H 9 (p)) = if Bq p^\qibb\B q (j))) = —if Bq 
the branching ratio of B S) d —> £ + £~ decays can be expressed as 


m 


B a 


m b + m q 


(3.32) 


B(B q 


0 + 


T Bq G%m 
87 r 5 


W£| [|Psm|-+ |V S m r] f 2 Bq m Bq m 2 t .l 


4 m 2 0 


m 


B, 


= J3(B q ^£+r)s M [\P\ 2 + \S\ 2 ], 


* [|P | 2 + |^| 2 ], 

(3.33) 


with 


D _ /"<SM 1 d'/ Bh b ^ySM q _ l-\ 

-rSM — '— ' l 0 + 0 2 ~ _ h/p , OSM — \ I t — 


277% m b + m Q 


4m| m b cSM 

m 2 Bq 2m^ m b + m q s 


and 


P = 


5 = 


1 


a/|-Psm | 2 + I^smI 2 
1 

v I-^SM| 2 + IS'smI 2 


^ , m B q rn b 

L 10 + o—2-1-°-F 

2 mp m h + m q 


= \P\ e ilpp , 


' 1 - 


4m < m ‘ Cs ee |S| . 


rrr Bi 2rnf v rn b + m, 


(3.34) 


(3.35) 

(3.36) 


While both Cf, M and Cf M are negligibly small, and hence Psm = Cf^, Ssm = 0, we keep their 
contributions in the numerical analysis both within the SM and in the MW model. In a more 
generic case, both P and S can carry nontrivial CP-violating phases p B and ps■ However, 
even in models with comparable Wilson coefficients, the contributions from 0$ and Op are 
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suppressed by a factor of /rn^y with respect to that from O w . Therefore, unless there were 
large enhancements for C$ and Cp, the coefficient Cio still provides the dominant contribution 
to the branching ratio. 

To compare with the experimental data, the effect of B q — B q oscillations should be taken 


into account, and the resulting averaged time-integrated branching ratio is given by 50,55 

1 + ^Ar Vq 


B(B q £+£-) = 


1 -y 2 a 


B(B q ^ £+£-), 


(3.37) 


where y q is related to the B q -meson decay width difference AT g , 

AT„ 


y q = 


r L- r tf 


n + n 


2 r„ 


(3.38) 


L'H ^ q 

with T q H{L) and T g = r^ 1 denoting the heavier (lighter)-eigenstate and the average decay widths, 
respectively. The time-dependent observable ^Ar can be, in the absence of new CP-violating 


contribution to the B q — B q mixing, expressed as 55 


^Ar = 


\P\ 2 cos 2 tp P — | S\ 2 cos 2 ip s 

\ P \ 2 + l ^| 2 


(3.39) 


In order to take into account the sizable AT S effect, it is convenient to introduce the ratio 55 

1 - Vs COS 2 (ps 


R = B{B s ^£+£ 


_ 1 + Vs cos 2 (pp | _ 12 j- — ys '-I" | c 12 

B(B S £+£-) sm ~ r PPy 2 1 1 1 - y 2 s 


|S| 2 - 


(3.40) 


It should be noted that both ^Ar and R can be directly extracted from the experimental 
observables [55]. Because of the negligible width difference in the B d system, however, the 
approximation B(B d —» £ + £~) ~ B(B d —* £ + £~) will be assumed throughout this paper. 


3.3 Z —y bb decay 

The Z-pole observable R®, defined as the ratio of the partial decay widths of the Z boson 
decaying into bottom quarks and into all quarks, 


, 0 _ r(z -> bb) _ r(z -> bb) 


R u h = 


T (Z —y hadrons) V(Z —y qq) ’ 

q^t 


(3.41) 


has been measured to high accuracy by the experiments at CERN e + e collider (LEP) and 
SLAC Linear Collider (SLC) [56], and plays an important role in precision tests of the SM, as 


well as in indirect searches for NP beyond it 57 . 
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The branching ratio R® is also quite sensitive to the NP models with an extended scalar 


sector due to the exchange of the additional charged and neutral scalars 22,58 . Specific to the 


MW model, the additional contribution to R® is dominated by the Z-penguin diagrams involving 
the charged colour-octet scalars, where the corresponding two-loop QCD corrections have been 


calculated for the first time by Degrassi and Slavich 22 . Following the same notations as in 


Ref. 22 , we have 


E ( gl ) + W 

-1 = i + 1 

Rh " 


7.1 \ 2 


Kn 


K h 


(3.42) 


{9l) +{9rY 

where g q , L are the left-handed and right-handed Zqq couplings, and the factor K q contains 
the QCD, QED and quark-mass corrections. The explicit expressions for these quantities can 


all be found in Ref. 22 


3.4 B — > X s 7 , B — > K* 7 and B —y p'j decays 

Both the inclusive B —> and the exclusive B —* V'f (V is a light vector meson) decays 

are induced by the quark-level b —» s(d )7 transitions, which contribute in a significant manner 


to current bounds on masses and interactions of possible charged scalars in NP models 59,60 


This is because the photonic-penguin diagrams with the charged-scalar exchanges contribute 
to these processes at the same level as those mediated by the SM W ± bosons. Besides the 
branching ratio B(B —y W/y), for which a good agreement has been observed between the 


experimental measurements |26,27 and the theoretical predictions 61 , we shall consider the 
two isospin asymmetries 


62,63 


r(R° k*° 7 ) - r(R+ ->• k*+ 7 ) 

f {B° -)■ K*° 7) + f (B+ K*+ 7 ) ’ 

(3.43) 

f {B+ p + 7 ) 

2 r(R° -»• p° 7 ) 

(3.44) 


A(py) = 


which exhibit a different dependence on the model parameters compared to B(B —> X s ~/) and 


could, therefore, provide complementary information on the model considered 121.64-661. 


I 11 the MW model, the additional contribution to b —» s(d )7 transitions comes only from the 
photonic-penguin diagrams exchanged by the charged colour-octet scalars. With the approx¬ 
imation m S) d —> 0, there appears no additional low-energy operator besides the SM ones, and 
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the NP effect is driven by new additive contributions to the Wilson coefficients at the matching 


scale fj,w, which have been calculated up to NLO by Degrassi and Slavich 22 . Following the 
same notations as in Ref. [22], we write the matching Wilson coefficients as 

a s (nw) 


Ci(fJ>w) = C-°\nw) + 5C\ u \nw) + 


( 0 ), 


47T 




(i)/ 


(3.45) 


67 


68 


where C\ k \nw) denote the SM contributions (k = 0,1) and could be found in Refs 
while 8C\ k \nw) denote the charged colour-octet scalar contributions. The explicit expressions 
for 5C[ k \^w) both at the LO and at the NLO could be found in Ref. 


22 


The evolution of 


the Wilson coefficients from the matching scale /i\v down to the low-energy scale fib, which is 
the same as in the SM, can be found in Refs. 167/68 j. The charged colour-octet scalar effects 


on the observables B(B —» X-y), A(py) and A(iF*y) have been studied in Ref. 21 , which will 
be followed in this paper. 


4 Numerical results and discussions 

With the theoretical framework at hand and the input parameters collected in Table [lj we now 
present and discuss our numerical results in this section. The stage II projection in the table 
refers to an epoch with 50 ftP 1 LHCb and 50 ab -1 Belle II data, which corresponds probably 


to the middle of 2020s at the earliest 35 . Estimates of future experimental and theoretical 


uncertainties are taken mainly from Refs. 33-35,71 • 


4.1 SM predictions and experimental data 

With the input parameters collected in Table [lj our predictions for the observables studied in 
this paper are listed in Table |2j where the current best measurements 26 27,31 together with 


the precision expected in stage II 133 35 are also given. The theoretical errors are obtained by 


varying each input parameter within the corresponding range and adding the individual error 
in quadrature. We assume that the central values of the future measurements in stage II remain 
the same as the current ones except for A (py). For the latter, it is observed that the current 
data disagrees with the SM expectation at more than 2cr level, and we assume therefore the 
central value of the future measurement to coincide with the SM expectation for A(py). 
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Tabic 1: Central values and uncertainties for the input parameters used in our analysis. The 

assumptions entering the stage II estimates are described in the text, and the entries “id” refer to the 
value in the same row as in the “2014” column. 


2014 


\v us \ {Ka) 

0.2258 ±0.0008 ±0.0012 

35 

69; 

\V cb \ x 10 3 

41.15 ±0.33 ±0.59 

35 

69 

\V ub \ x 10 3 

3.75 ±0.14 ±0.26 

35 

69 

7 

(68.0±|;°)° 

35 

69 

f B . [GeV] 

0.228 ± 0.008 

49 

70 ; 

fs d [GeV] 

0.189 ±0.008 

49 

70 

fs 3 /fB d 

1.206 ± 0.024 

49 

70 ; 


m c 

1.275 ± 0.025 

26 

rh b 

4.15 ±0.03 

26 

pole 

m 

173.21 ±0.51 ±0.71 

26 

m h o 

125.7 ±0.4 

26 

sin 2 6 W 

0.23126 ± 0.00005 

26 

a s {m z ) 

0.1185 ±0.0006 

26 


Stage II 

0.22494 ± 0.0006 
42.3 ±0.3 
3.56 ±0.08 
(67.1 ± 1)° 
0.232 ±0.001 
0.193 ±0.001 
1.205 ±0.005 
1.286 ±0.010 
id 
id 
id 
id 
id 


35 


34 


34 


33 


34 


35 


35 


35 


71 


71 


71 



Taking into account the theoretical and experimental uncertainties, one can see that the SM 
predictions presented here agree with the corresponding experimental data within 2a error bars, 
and hence strong constraints on the model parameters are expected from these observables. 


4.2 Procedure in the numerical analysis 

The flavour observables discussed in this paper involve only four parameters of the MW model: 
rj ul rj d ,m S A and Ai- As discussed in Sec. |2j the parameter Ai is real due to the Hermiticity of 


the scalar potential 17 . For the Yukawa coupling parameters r) u and rid, we shall consider the 
following two cases: 

(i) both of Tj u and rid are real, which is referred as “real couplings” throughout this paper. 
Unlike in the case of the usual 2HDMs [7|, these two parameters are assumed to be 
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Tabic 2: SM predictions and experimental measurements for the observables discussed in this paper. 
The entries “id” refer to the value in the same row as in the “2014” column. 


Observable 

Exp. 

data 


SM prediction 


2014 


Stage II 


2014 

Stage II 

Am B o [ps _1 ] 

17.761 ±0.022 

26 

id 


17.568±l;i" 

I 7 S89 +0 ' 431 

I I .ooy_ 0 _ 429 

Am B o [ps" 1 ] 

0.510 ±0.003 

26 

id 


D KQQ+0-086 

U.OOO_q 085 

0 51 5+ 0 - 016 
u.oio_ 0 016 

B(B S — > /T+/T') [10“ 9 ] 

9 8+ 0 - 7 
^• 0 - 0.6 

31 

2.8 ±0.3 

33 


o 46^0-33 
O '^ u -0.30 

q 7 q + 016 
°- , y - 0.22 

B(B d — >■ fx+fi-) [10- 10 ] 

3.9±i;| 

31 

3.9 ± 1.4 

33 


n qq+0-16 
u - aa -0.16 

1 07 + °' 05 
i - u ‘-0.06 

B(B — >■ X sl ) [10- 4 ] 

3.43 ±0.22 

27 

3.43 ±0.21 

34 


3.i4i?:i? 

oil +0.25 
°- i± -0.26 

A(iC* 7 ) [10“ 2 ] 

5.2 ±2.6 

27 

5.2 ±0.3 

34 


4.15± 2 ;^ 

4.20+ 2 ;^ 

A(P7) [10- 2 ] 

-46 ± 17 

27 

-7.4 ± 1.7 

34 


-8.29 • ^ 

-7.38+l;|I 

R° b [10- 3 ] 

216.29 ±0.66 

26 

id 

215.86±8:8i 

id 


independent from each other. 


(ii) both r] u and rjd are complex, which is referred as “complex couplings” in the following. In 
this case, the independent Yukawa coupling parameters can be chosen as the magnitudes 
\rj v \ and r ] d |, and the relative phase 6 defined through r]*j] d = 

Constraints on the model parameters have been analyzed both theoretically by requiring the 


perturbative unitarity and the vacuum stability 20 , as well as phenomenologically by satisfying 


the current flavour, Higgs and EW precision data 17,21-25 . It is important to mention that 


the current experimental searches based on the dijet events by the ATLAS 72 and CMS 73 


collaborations exclude already the existence of colour-octet scalars with mass below ~ 2 TeV, 


while the four-jet searches at the LHC do not observe any signal at a much lower mass region 74 


75 . However, if the colour-octet scalar decays into more than two light quarks, top-quark and 
jet, or tt, the existing bounds on the colour-octet scalar mass would be different and even 


masses of order 400 GeV could not be excluded, as noticed in 17,24 . In addition, all these 


mass bounds depend on some assumptions about the underlying NP models. In this paper, 
therefore, we shall specify the model parameters within the following ranges: 


Ai G 


[-40,40], m S A e [300,1000] GeV, r lu e [-3,3], % e [-100,100] (4.1) 
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in the case of real couplings, and 


Xi G [-40,40], m S A G [300,1000] GeV, \r) u \ G [0,3], \r) d \ G [0,100], 6 G [-180°, 180°] (4.2) 


in the case of complex couplings. Then, we make an extensive random scan for these model 
parameters over the above constrained ranges to generate the initial samples. 

For the numerical analysis, we adopt the following procedure: for each sample generated 


from the parameter ranges specified by Eqs. 4.1 and 4.2, we give the theoretical prediction 
for an observable together with the corresponding theoretical uncertainty. If the obtained 
theoretical range for the observable has overlap with the 2cr range of the experimental data, 
the sample is regarded to be allowed. Following this procedure, we get the final survived ranges 
for the model parameters. To incorporate the theoretical uncertainty, we use the statistical 


treatment based on frequentist statistics and Rht scheme 76 , which has been implemented in 


the CKMhtter package 69 . Here the main observation is that, while the experimental data 


yields approximative^ a Gaussian distribution, the theoretical calculation for an observable 
does not. The latter depends on a set of input parameters like form factors, decay constants 
and Gegenbauer moments etc., for which no probability distribution is known. The Rht scheme 
assumes no particular distribution for the theory parameters, only that they are constrained 
to certain allowed ranges with an equal weighting, irrespective of how close they are from the 
edges of the allowed range. In addition, for simplicity, the relative theoretical uncertainty is 
assumed to be constant at each point in the parameter space. This is a reasonable assumption, 
since the main theoretical uncertainties are due to the hadronic input parameters, common 
to both the SM and the NP contributions. Therefore, the theoretical range for an observable 
at each point is obtained by varying each input parameter within its respective allowed range 
and then adding the individual uncertainty in quadrature. This procedure has already been 
adopted in Refs. 121164,771. 

In order to investigate the capability of the MW model in explaining the flavour physics 
data listed in Table [2] and to check the validity of the samples satisfying all the experimental 
constraints, we construct the y 2 function following the method proposed in Ref. 

(or p - of w ) 2 


78 


x 2 = 


E 

zG(o6s.) 


K exp ) 2 + (^ MW ) 2 ’ 


(4.3) 
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where 0® xp and 0^ w denote the central values of the experimental measurements and the 
theoretical predictions for an observable *, with a“ p and a^ IW being the corresponding experi¬ 
mental and theoretical errors, respectively. Using the method, one can obtain the y 2 value for 
each point in the parameter space, and then the minimal y 2 value. 


4.3 -B s ,d — B s d mixings within the MW model 


As mentioned already in Sec. 3.1, the B sd — B s>d mixings within the MW model involve only 
the four Wilson coefficients Cy RL , Cy[ L , C RRR1 and C RRR2 . Taking m S A = 500 GeV as a 
benchmark, we obtain numerically: 


CyLL (fh>) 
(») 


= 0.05|t/ u | 2 + 0.02|r,„| 4 + ( — ) • 1(T 3 1.02< % + 2.35|rj„| 2 + 0.22|i,„| 4 , (4.4) 

m w 1 


Csrr, 1 (/ b >) _ 

A/LL (/h>) V m W 


mb ' .10- 2 


1.05 rf u r] d - Q.521^1 2 + O.U( V * uVd ) 2 - O.U V * uVd \ri u \ 2 + 0.03|^| 4 


(4.5) 


CsRR ,2 (A(&) _ 
CyLL M 


mb ' .10- 4 


3.22|^| 2 - 6.55 rj* u rj d - l.85( V * uVd ) 2 + 1.85^1^1 2 - 0.36|*/ n | 4 


(4.6) 


for the short-distance Wilson coefficients, and 


(. B s \n[g Bl ~ 2 \B s ) = <j 5.77 - 0.22* + (0.29 - 0.01*) \r]u{ Z + 10” 2 • (8.77 - 0.34z) \q. 


m b 

m w 


• 10 


-2 


(4.27 - 0.17*) rf u r ld + (3.75 - 0.15*) \ Vv 


- (0.51 - 0.02*) (■ rf u rj d f + (0.51 - 0.02*) vMuf 


+ 0.03KI 4 f x 10 -12 GeV, 


(4.7) 


for the off-diagonal matrix element. From these results, we make the following observations: 


• The theoretical expressions for the short-distance Wilson coefficients are independent of 
the parameter Ai, and hence the mixing observable Am^o puts no constraints on the 
trilinear scalar coupling g S A S B h0 . 
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• The Wilson coefficients Csrr \ an d ^srr 2 are suppressed by a factor of rn^/rn^ with 
respect to Cy[ L . So the dominant NP contribution to the mass difference Am B o comes 
from the operator Q\ LL and is always constructive to the SM one. 

• Except for the case with much smaller rj u , the terms containing 77*77,2 in the off-diagonal 
matrix element are suppressed by a factor of ml/m^y. So the constraint on the parameter 
rid from the mixing observable Am B o alone is quite weak. 

• With the 77*77,2 terms neglected, the off-diagonal matrix element and hence the mass differ¬ 
ence Am B o depend on the parameter r} u only via \r) u \. So the case with complex couplings 
can be inferred trivially from that with real couplings. 

The bounds on the model parameters derived from Am B o are shown in Fig. [4] in the case 
of complex couplings. As expected, the regions with large \rj u \ are already excluded due to 
the good agreement between the SM prediction and the experimental measurement for Am B o. 
The current bounds on the model parameters can be significantly improved by the stage II 
projection. There are, however, almost no constraints from Am B o on the coupling rj ( i and the 
phase 9. Except for being a little bit weaker, the bounds on the model parameters derived from 
Am B o are similar to that from Am B o, and will not be shown here. 


4.4 B s d —V n + n within the MW model 


The B s —> /i + fi decay has been observed with a statistical significance exceeding 6cr and, 
furthermore, a 3<v evidence for the Bd —> /i + /i~ decay is obtained by the CMS and LHCb 


experiments 31 . Both measurements are statistically compatible with the SM predictions and 


place, therefore, stringent constraints on physics beyond the SM 28,50 . We now investigate 
the constraints on the model parameters from the branching ratios of these two processes. 


As mentioned already in Sec. |3.2[ the rare leptonic B-meson decays are dominated by the 
Wilson coefficient Cio both within the SM and in the MW model. To improve the accuracy of 


the SM prediction, we use the fitting formula for Cf^ 1 given by Eq. (4) of Ref. 32 , which has 
been transformed to our convention for the effective weak Hamiltonian 1281 


^-ySM _ 
L '10 — 


-0.9604 


M t 

1.52 

&s(Mz) 

-0.09 

+ 0.0224 

[ M t 

0.89 

Oi s (M z ) 

173.1 GeV 


_ 0.1184 

173.1 GeV 


0.1184 _ 


-0.09 


(4,8) 
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Figure 4: Bounds on the model parameters from Arn^o in the case of complex couplings. 

As the NP contribution to C'io in the MW model involves only the parameter rj u and the mass 
m S A, stringent bounds on them are expected from these processes. For a benchmark value 
m S A = 500 GeV, the coefficients P and S read numerically 

.2 


777/ / 

P = -1.003 - 0.138ffiJ 2 + —(0.586 + 0.003W + 0.068ffi, 

m z w V 


S = 


m 


2 r 

b s 


m 


W L 


-0.306 - AG 0.020^ - 0.006|^| 2 ) - 0.075^ + 0.008|7/, 


from which, we make the following two observations: 


(4.9) 

(4.10) 


In Eq. (4.9), the term proportional to |?/ n | 2 corresponds to the NP contribution to G 10 , 
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Figure 5: Constraints on the model parameters from B(B S —?• ) in the case of real couplings. 

and is not suppressed by the factor m 2 Rn jm‘\ ¥ . It contributes always constructively to 
C',q T , irrespectively of whether r] u is real or complex. 

• The terms involving r]d and Ai are all suppressed by a factor of rri^Jrri^ and their effects 
are, therefore, negligibly small unless when r]d or Ai is abnormally large. 

Under the constraints from the branching ratio B(B S —> /i + /i ~), the allowed parameter 
spaces are shown in Figs. [5] and [6j corresponding to the case of real and complex couplings, 
respectively. As expected, the observable can provide strong bounds on \rj u \ and m S A, but not 
on r/d and Ai, irrespective of whether they are real or complex. In the stage II epoch, as both 
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Figure 6: Constraints on the model parameters from B(B S —> ji ) in the case of complex couplings. 


the experimental and the theoretical uncertainties will become much smaller and, furthermore, 
due to the constructive interference between the SM and the NP contributions, the parameter 
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Figure 7: Constraints on the model parameters from B(Bd —> g + g ) in the case of real (the upper 
panel) and complex ( the lower panel) couplings. 

spaces allowed by the observable B(B S —y g + g~) will be reduced significantly. 

As shown in Fig. [tJ, the constraints from B(Bd —» p + /i~) are much weaker compared to that 
from B(B S —>■ /i + /z~), which is understandable due to the following two facts: 

• Besides the much larger error bars, the experimental data on B(Bd —> /i + p~) is almost 
2a above the corresponding SM prediction. An opposite behaviour is, however, observed 
for B(B S —> 

• For both of these two processes, due to the constructive interference between the SM and 
NP contributions, the NP effect can only enhance the theoretical prediction. 


4.5 R® within the MW model 


For the Z —» bb decay, the NP effect resides only in the left-handed and right-handed Zbb 
couplings g b L R and, in the assumption of negligible oblique corrections due to the second scalar 
doublet, only the parameters rj u , r) ( i and m S A are involved 22 . Taking m S A = 500 GeV as a 
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Figure 8: Constraints on the model parameters from R® in the case of real (the upper panel) and 
complex (the lower panel) couplings. 


benchmark, we get numerically 

G g b L ) 2 = 0.18 - 5.5 x 10- 4 |?/J 2 + 2.34 x Kr 7 |^| 4 - 5.5 x 1(T 10 M 2 - 2.9 x 10~ 10 |7 7 ^ 7rf | 2 , (4.11) 


(g b R ) 2 = 0.006 - 1.9 x 10~ 7 \r] u \ 2 - 3.1 x 1Q- 8 |^| 4 - 2.6 x 10- 8 |%| 2 + 1.9 x lO- 11 ^ 2 . (4.12) 


Together with Eq. (3.42), the following features are observed: 


• The ratio R® is sensitive to the parameter g u , but not to rjd- This is due to the fact that 
the terms controlled by rjd are suppressed by a factor of m 2 /m 2 with respect to the ones 
controlled by rj u . 

• The NP effect is dominated by the \rj u | 2 terms and contributes destructively to the SM 
one. So the contribution from the MW model will always decrease the value of R®. 

• As both of the couplings g b L and g b R depend only on the magnitudes of g u and r)d, no 
information on the phase 9 can be inferred from the ratio R®. 


In Fig. [8j we present the constraints on the model parameters from R®. As expected, there 
are almost no constraints on 7]d, irrespective of whether it is real or complex. In both cases, 
only small values for \rj u \ are allowed by the current data. It is also noted that the constraints 
in stage II remain almost the same as that in 2014 and are, therefore, not shown here. 
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Figure 9: Constraints on the model parameters from B —> Ms 7 in the case of real couplings. 


4.6 B{B —y Xg'y), A(K*~/) and within the MW model 

The dominant contributions to the radiative b —» s(d )7 transitions in the MW model are 

related to the effective Wilson coefficients of dipole operators and Cf® (//*,). In the 

leading logarithmic approximation and taking m S A = 500 GeV, we get numerically 

= 1 - 0.11<% + 0.012|%| 2 , (4.13) 

CfM/Cf^) = 1 + 0-1377:% - 0.022|%| 2 . (4.14) 


It is observed that the dominant NP contribution comes from the term proportional to the 
combination r\* u r\d when r\ u and % are comparable. In the LO approximation, the branching ratio 
B(B — y Xg'y) is known to be proportional to \Cf^ (//.&)| 2 . This observable can, therefore, provide 
stringent constraints on the combination rj^rjci- This is shown in Figs. [9] and 10, corresponding 
to the case of real and complex couplings, respectively. Based on these plots, we make the 
following observations: 

• As shown in Fig. [9j in the case of real couplings, there are two allowed regions under the 
constraint from B(B —» A s y). The region close to the axes corresponds to the case when 
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Figure 10: Constraints on the model parameters from B —> X s j in the case of complex couplings. 

the NP contribution is small and constructive to the SM one. The other region, in which 
large and same-sign values for r/ u and 7] d are allowed simultaneously, corresponds to the 
case when the NP contribution is destructive to the SM one and makes the coefficient 
Cy ff (/i 6 ) sign-flipped. It is also noted that, the regions with simultaneously large values 
for rj u and rj d but with opposite signs are already excluded. 

• As shown in Fig. [To^a-b) , the current data on B(B —>• A" s y) also gives strong constraints 
on the parameters \rj u \ and \r)d\ in the case of complex couplings. Furthermore, the 
interference between the SM and NP contributions depends on the phase 9, as shown in 
Fig. [To|e-f). When 6 « ±180°, the interference is constructive and only a small region 
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Figure 11: Constraints on the model parameters from A(/\* 7 ): (a-b) for real couplings, while (c-e) 
for complex couplings. 


with smaller 1 77 * 77 , 2 1 remains. While for 6 « ±0°, the interference becomes destructive 
and there exist two allowed regions, corresponding to the case with relatively small NP 
influence (the lower region) and the case when the NP contribution is about twice the 
size of the SM one (the upper region), respectively. 

• As shown in Figs. |9](c-d) and[!o|(c-d), the allowed values for the combination r] u 7j d (in the 
case of real couplings) and 77 * 77,2 | (in the case of complex couplings) increase with the 
mass of the charged colour-octet scalars, as they should be. 


For the exclusive B —> V'y decays, the branching ratios could not provide further con¬ 
straints on the model parameters compared to that of the inclusive B —> Ah, 7 decay. However, 


the two isospin asymmetries defined by Eqs. (3.43) and (3.44), which show a different depen¬ 
dence on the model parameters from that of the branching ratios, may provide complementary 


constraints 1211. Under the constraints from A(A"*y) and A(p 7 ), we show in Figs. 11 and 


12 the allowed regions for the model parameters. From these plots, we make the following 


observations: 


• As for A(/F* 7 ), the 2014 data has almost no constraints on the model parameters and 
are, therefore, not shown here. In the stage II projection, however, because of the re- 
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Figure 12: Constraints on the model parameters from A(p'y): (a-c) for real couplings, while (d-f) for 
complex couplings. 


duced uncertainties in both experimental measurements and theoretical predictions, the 
constraints from this observable become significant, which results in an upper bound on 
the combination | r[*r ] d |, as shown in Figs. [IT|b) and [fijd). Strong correlations are also 
observed in the q u — r} d , m S A — r} u rj d (in the case of real couplings), as well as in the 
\v u \ ~ \Vd\, m s a — \r)*rjd\ and 6 — \ r}^r} d \ (in the case of complex couplings) planes, as shown 
in Figs. [lT|(a), [TT|(b) , [IT|(c) , [lT|(d) and [IT|(e) , respectively. 


As for A (py), the 2014 data provides already an upper bound on \rj u rjd\ in the case of 


real couplings, as shown in Fig. 12 b). Moreover, under the constraint of the 2014 data, 
strong correlations exist in the r] u — r] d , m S A — r] u r]d (in the case of real couplings), as well 
as in the m S A — and 0 — \rj*rjd\ (in the case of complex couplings) planes, which are 

plotted in Figs. |l2|(a), 12 'b),[l2|(d) and[l2|'e), respectively. In the stage II epoch, however, 
only weak bounds on the model parameters can be obtained, as shown in Figs. [I2|(c) and 


12 f). This is because the allowed range for NP contribution is much larger than that in 


2014, when the central value of the future measurement is chosen to coincide with the SM 
prediction and the interference between the SM and NP contributions is constructive. 


• As can be seen from Figs. |To|e-f), Fig. [Tlje) and Fig. |l2|e-f), the branching ratio B(B —>■ 
X s y) and the two isospin asymmetries A(A'*y) and A (py) exhibit a different dependence 
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Figure 13: Combined constraints on the model parameters in the case of real couplings, 
on the phase 9. So a combined constraint from these observables should be more stringent, 


which will be discussed later 21 . 


4.7 Combined analysis within the MW model 


Combining all the constraints from the observables discussed in the previous subsections, the 


final survived parameter spaces of the MW model are shown in Fig. 13 and 14, corresponding 
to the case of real and complex couplings, respectively. From these plots, we make the following 
observations: 


I). In the case of real couplings 

• The parameter r] u is strongly bounded by Am B oB(B Sj d —> /J, + n~) and with 
large values for it being already excluded. The most stringent constraint on rj u comes 
from the branching ratio B(B S —> /i + /i _ ). 
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All the observables considered put no bounds on the individual parameter r) d . How¬ 
ever, the combination r} u rj d is significantly bounded by B(B —>■ A s y), A(A'*y) and 
A(py), with only small values for it being allowed, as shown in Figs. [l3|^d) an d mi(g)- 
Strong constraints on the combination r] u ‘rj d can also be seen in the correlation be¬ 
tween rj u and r ] d , as shown in Figs. [I3fa) and 13 'b). 


There is no limit on the individual parameter Ai, which can be seen from Figs. |l3[e) 
h). This is mainly because the parameter Ai is relevant only to B(B Sjd —)■ 


and 


13 


p + p ) and the terms containing Ai are, however, suppressed by a factor of m% d /rn^ v 
as discussed in Sec. 14.41 

Bounds on the mass of the charged colour-octet scalars are always accompanied by 
the Yukawa coupling parameters t] Ul r] d or their combination, which results in strong 
correlations between rj u and m S A or between rj u r] d and m S A, as shown in Figs. 13Jc-d) 


and p~3| f-g). It is also noted that the observables B(B —>• A s y), A(A'*y) and A(py) 
provide strong correlations only in the r] u r] d — m S A plane, while the observables 
Am B o , B(B S)d —> /i + /i~) and A° have significant effects only on the ii u — m S A plane. 


II). In the case of complex couplings 


• Similar to the case of real couplings, the magnitude Tj u is strongly bounded by 
Amgo^, B(B S:d —> /U + /i”) and A°, with large values for it being already excluded. 
The most stringent constraint on \rj u \ comes also from B(B S —> 


• No bounds on the parameter \r) d \ are provided from the observables discussed in this 
paper. However, the combination is significantly bounded by B(B —>■ AAy), 

A(A'*y) and A(py), as shown in Figs. [I4|(e-f) and[l4|(h-i). The strong constraints on 
|y* 7 d| from these observables can also be seen in the correlations between \r] u \ and 
\rj d \ shown in Figs. [l4](a) and[l4^c). 


As can be seen from Figs. 14 g) and 14), there is no limit on the parameter Ai, 
which is due to the same reason as in the case of real couplings. It implies also that 
there is no correlation between \r] u \ and Ai under the combined constraints. 
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Figure 14: Combined constraints on the model parameters in the case of complex couplings. 


There exist strong correlations between \r) u \ and m S A or between \r)*rjd\ and m S A, 


which can be seen from Figs. 14 [b), |l4[d-e) and 14(h). Here complementary con¬ 
straints on the model parameters from the observables B(B —> X s 7 ), and 

A(py) are crucial. 

As the phase 9 is always associated with the combination \r)*/ld\, strong correlations 


between them are observed in the 6 — |?/*?^| plane, which can be seen from Figs. 14 f) 
and[l4|(i). 


We keep the initial samples generated within the ranges specified by Eqs. (4.1) and (4.2), 


and calculate the corresponding y 2 associated with each of the samples. It is found that, for 
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Figure 15: The Ay 2 distribution of the data samples allowed by the experimental data on the observ¬ 
ables discussed in the paper. The different colours correspond to different values for the cumulative 
distribution function. 


the 2014 data, y 2 /d.o.f. is 10 within the SM, while the minimal value for it equals 9.5 (for real 
couplings) and 6.8 (for complex couplings) in the MW model. Here d.o.f. denotes the number 
of free parameters in a considered model and, specific to the MW model, it is 4 in the case of 
real and 5 in the case of complex couplings, respectively. It is noted that Xmin f° r d ie case °f 
complex couplings is much smaller than that for the case of real couplings. This is mainly due 
to the different values of d.o.f. in the two cases. It is, therefore, concluded that the MW model 
is well suited to explain the current data on the observables discussed in the paper, with the 
resulting y 2 being significantly smaller than the SM counterpart. 

To explore the validity of the final data samples obtained by imposing all the constraints in 


the way mentioned in Sec. 4.2, we calculate the Ay 2 associated with each point in the samples 


and show in Fig. 15 the Ay 2 distribution of the points. It is concluded that almost all the 
points lie in the 99.7% confidence intervals, indicating that the procedure adopted in this paper 
is reliable and conservative. 


4.8 Correlations between different observables in the MW model 

Up to now, only B(B s d —* /i + /r _ ) have been measured among the observables for B sd —>■ 
/i + /r“ decays. In order to gain further insights into the model parameters from the future 
measurements, we present in Fig. [16] the correlations between the observables dAr and R 
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Figure 16: Correlation plots in the R — ^Ar plane under the combined constraints. The allowed 
regions are shown by green (dark grey) and grey (light grey) points, which are obtained in the case of 
real and complex couplings, respectively. 


defined respectively by Eqs. (3.39) and (3.40), in the survived parameter spaces under the 


constraints discussed in previous subsections. 

As some of the observables discussed in this paper still suffer large uncertainties, we also 
investigate the correlations among them under the constraints from the other observables, which 


is shown in Fig. 17 For simplicity, we do not consider the theoretical uncertainty at each point 


in the parameter space. Given the assumption that the relative theoretical uncertainties of these 
observables are constant at each point in the parameters space, we show the SM predictions 
with the corresponding theoretical range in these plots, which can be applied to each point to 
account for the theoretical uncertainty at that point. 


From the correlation plots shown in Figs. 16 and 17, we make the following observations: 


In the MW model, large deviations from the SM prediction for the ratio R are still 
allowed. The observable ^Ar remains, however, almost independent of the NP effect, 
because the NP contribution to ^Ar is suppressed by a factor of In the stage II 

projection, the allowed range for R becomes much smaller due to the stringent constraint 
from B(B S —> /i + p~), which has already been discussed in Sec. 


4.4 


It is also noted that 


the NP effect always increases the value of R with respect to the SM prediction. 

Under the constraints from Am B o d , B(B —> Ay/), A(A'*y) and A(py), the correlation 

between B(B S —> /i + /i _ ) and B(B d —> /i + /i”) is very strong, as shown in Figs. 


17 


and mid). The behaviour of the correlation is independent of whether the parameters 
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Figure 17: Correlation plots among B(B Sd -> /r + /u~), A(A'*y) and A(py). The dashed lines denote 
the experimental data with 2cr error bars, while the SM predictions with the corresponding 2cr range 
are shown by the blue (dark) cross. This cross is also applied to each point to account for the theoretical 


uncertainty at that point. The other captions are the same as in Fig. 16 


are real or complex, because both of these two observables are dominated by the terms 
proportional to r ] u | 2 , as discussed already in Sec. 

• Under the constraints from Am B o ^ B(B —> X s y), R® and B(B Std —> /i + /r~), the observ¬ 
able A(A^y) is strongly correlated with A(py), as shown in Figs. [I7|(b) and [lT](e) . This 
is mainly because both of them are sensitive to the value of ^uVd (for real couplings), and 
|?7*?7d| and 9 (for complex couplings). The correlation between A(A'*y) and A(py) in the 
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case of real couplings is, however, different from that for complex couplings. 


As discussed already in Sec. 4.4, the observable B{B S —» /z + /i ) in the stage II projection 


can provide stringent bounds on Tj u |, restricting the theoretical prediction for B(B d —> 
to a very narrow range under the constraint from B(B S —>- /j, + /j,~ ), which can be 
clearly seen from Figs. [l7^c) , 17 f), [f^h) and 17'k). 


• As the observables A(A'*y), A(py) and B(B S)d — > /x + /r~) exhibit a quite different de¬ 
pendence on the model parameters, the parameter space allowed by A(AT*y) or A (py) is 
complementary to that allowed by B(B sd —> p + p~). This fact results in mild correlations 
between either of A(AT*y) and A(py) and either of B(B S —> /i + p _ ) and B(B d —> p + p~), 
as shown in Figs. |l7|^c), Qf) and[l7[g-l). 


4.9 Comparison with collider constraints within the MW model 

The phenomenological aspects of the MW model at hadron colliders, including the single and 
pair productions of these colour-octet scalars, their decays, as well as their implications in Higgs 
data have been studies extensively in the literature [17,24,25]. It should be noted that the LHC 
signatures of the colour-octet scalars depend strongly on the model parameters, such as the 
masses of the colour-octet scalars, as well as the values of rj u , j] d and A*. As noticed already in 


Refs. 17,24 , if the colour-octet scalar decays into more than two light quarks, top-quark and 
jet, or tt, the existing bounds on the colour-octet scalar mass from ATLAS and CMS would be 
different and even masses of order 400 GeV could not be excluded. 

The low-energy processes considered by us, on the other hand, involve only the parameters 
Vu, Vd, an d Ai. It is found that these low-energy flavour observables could put stronger 
bounds on the parameter and the combination rj u rj d , especially with the observables achiev¬ 
able with 50 fb -1 of LHCb and 50 air 1 of Belle-II data. Furthermore, under these indirect 
constraints, strong correlations are observed between t] u and m S A and between r] u ri d and m S A. 
All these indirect constraints can provide very important and complementary information on 
the model parameters, which is useful for the direct searches of the colour-octet scalars at 
hadron colliders. 
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5 Conclusion 


In this paper, we have performed a complete one-loop computation of the short-distance Wilson 
coefficients for B S)d — B sd mixings and B s>d —> £ + £~ decays within the MW model. It is found 
that, in order to get a gauge-independent result, the external momenta of the heavy quarks 
inside the mesons should be taken into account, and the heavy-quark masses should be kept 
up to the second order. 

Based on these calculations, combined constraints on the model parameters have been de¬ 
rived from the current flavour data, including the B s d — B s d mixings, B s d —> B —> Ah-y, 

B —> K* 7 , B — > py, and Z —>• bb decays. The future sensitivity to the model has also been 
explored in the observables achievable with 50 fb -1 of LHCb and 50 ab _1 of Bcllc-II data. Our 
main conclusions are summarized as follows: 

• The flavour observables exhibit a different dependence on the model parameters and 

could, therefore, provide complementary bounds on the parameter space. The observables 
Am B oB(B Std —> and R® are shown to be only sensitive to the parameter rj u . 

Furthermore, there exist strong correlations between rj u and m S A or between \rj u \ and 
m S A. On the other hand, the observables B(B —> X s j) , A(/l* 7 ) and A(py) are found to 
be only sensitive to the combination rf v rj d . Strong correlations between rj u rj d and m S A or 
between rj*rj d and either of rn S A and 6 are also observed. 

• The correlation between the observables ^Ar and R in the B s —* /i + fx~ decay, which is 
characterized by the sizable decay width difference of B s system, is investigated within the 
survived parameter space from all the constraints. While large deviations from the SM 
prediction are still allowed for the ratio R, the observable ^Ar remains almost independent 
of the NP contribution. It is also noted that the NP effect always increases the ratio R 
with respect to the SM prediction. 

• We have calculated the \ 2 value associated with each of the generated samples and con¬ 
cluded that the MW model in the case of complex couplings is more suitable to explain 
the current flavour data discussed in this paper than in the case of real couplings. It 
is also concluded that the MW model in both cases is capable of explaining the current 
data, with the resulting Xmin being significantly smaller than that of the SM. 


38 



• Since some of the observables still possess large uncertainties, we have also investigated 
correlations among the observables A(K*^y), A(py) and B{B sd —> within the al¬ 

lowed parameter space. Some of the correlations are strong and could provide further 
insights into the MW model, once more precise experimental measurements and theoret¬ 
ical predictions are available in the future. 

In summary, the results obtained in this paper are helpful to deepen our understanding 
of the MW model, and could be complementary to the analysis performed at the high-energy 
frontier of the LHC. With the promising progress expected from the LHCb and Bellc-II, and 
the improved theoretical predictions, more detailed information on the model can be inferred 
from these low-energy processes. 
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A Relevant coefficients for B s .d — B s ^ mixings 
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Here we present the explicit expressions for C VL ' L , C VL ’ L , C VLL , C SRR1 , C SRRXl C SRRAl 


C SRR% C SRR ,2 and C SRR% appearing in Eq. Q: 
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6 (xs - 1) (xt - 1) (x c (xs + 1) - xs) (xs - Xt) ’ 


6(x t - 1) (xt - x c ) (xs - xt)" 


(A.16) 


C 


NP,ct 


19x fe x c x|xt (x s - 3x t ) 19x fe x c xsXj (x s - 3x t ) 

-j In x c H- t In x+ 

\4 ° mo / o _ __ \ .. \4 1 


SRR\\r)u I 208 (4x c — £5) (X5 — X t ) 


108 (3x c - Xt) (xs - Xt) 
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19 x b x c x 2 t (—22 xsx t + 5x| + 5xf) 
648(x s - xtf ( x c (2x s + 3 x t ) - x s x t ) ’ 


(A.17) 


£jNP,ct 


'skrx,\ vu\ 2 = ^ 32 ^ c ’ XS ’ Xt ) 111 XS ~ ^ 33 ^ c ’ XS ’ X *) ^ 34 ( Xc ’ X *) 111 Xt 


+ /35 (®c, as, a*, Xt) / 36 (x c , x s , x b , x t ) + f 37 (x c , x s , x b , Xt) 


(A.18) 


c 


NP,tt 

SRRl,ridri* |r/ u | 2 


19x fe xf (x s + xt) 19x fe x 2 (x s + xt) 19x 6 x 2 

-o— In xq H-o— m x t H- 

\3 ° i ~ \3 1 


18(x s - x t y 


i8(x s - x t y 


9(x s - xt) 


2 » 


(A.19) 


C 


NP,tt 


SRRl,(r) d ri*) 18 (x^ — X t ) 


19x b x 2 (xs + Xt) 19x 6 x 2 (xs + x t ) . 19x fe xf 

- : -5— m xc- 1 -5— m Xf- 

3 ° 1 o/_. _ \3 z 


18(x s - x t y 


9(x s - Xt) 


2 > 


(A.20) 


s~^N P,tt 
^ SRRl,>ldVt 


5x b x s x 2 (xs (xt + 1) - 2Xt) , 5x fe xf (-2x 5 x t + x| - (x t - 2) xf) 

- 3 - In Xs H-—--7s- In x t 


3(x s - l) 2 (x s - xt) 3 3(x t - l) 3 (x s - xt) 3 

5x fe xf (x| (xt - 3) + xs (—3xf + 6x t + 1) + (x t - 3) x t ) 

6 (x s - 1) (x t - l) 2 (x s - xt) 2 


(A.21) 


/~^N P,tt 
U SRR1, \ Vu \ 4 


19x;,xf (x 3 + xf — 3x|x t — 3x^xf) ^ 19xbxf (x| + xf — 3x|x t — 3xsx 2 ) ^ 


108(xs - xt) 5 

19x fe xf (—22x§xt + 5x| + 5xf) 
324(x s - xt) 4 


108(xs - xt)^ 


(A.22) 


^SSi.i^i 2 = ^ 3s ^ c ’ a:s ’ Xt ) ln + ^ 39 xs ’ ln Xt 


- ho (x c , x s , X b , Xt) / 4 1 {x c ,x s ,x b ,x t ) + / 42 (x c ,x s ,x b ,x t ) + / 43 (x c ,x s ,x b ,x t ) , (A.23) 


C 


NP.cc 


’SRRl,r} d v*\vu\ 2 ^SRRl,(r, dV *) 2 U ’ ^SRRl, Vd ,^ U > L SRRl,\ Vu \ 4 U ’ ^SRRl,\i lu \ 2 




NP,ct _ 3 (~fNP,ct siNP,ct _ 21 ^yNP,ct iNP,ct. 

U SRR2,r] d riZ ~ _ 20 ^ SRRltfdVV ° SRR2,\ Vu \ 4 ~ ~ 75 L ' SRRl,\ri u \ 4 ’ °Si?R 2,|^ u | 2 


^NP,ct _ JA r NP,ct ^NP,ct _ 21 „ATP ;C t 

°SRR2,^r;S|r , u | 2 ~~ 76°Si?Rl, % <|^| 2 ’ L 'SRR2,{ Vdn *) 2 ~ SRRl,( VdV *) 2 ^ 


6 ,-iNP,ct 
' 20 U SRRl,\ Vu \ 2 ’ 

(A.25) 


s~iNP,tt 
^SRR2,rj d ri* 


3 x-y NP,tt y^l 

20 SRRl ^dvV l 'srr 2,\ Vu \ 4 


iNP,tt 


21 ^,NP,t.t 
7g U SfiRl,|r; u | 4 ’ 


s~-iN P,tt 
U SRR2,\ nu \ 2 


6 s~iNP,tt 
'20 SRRl,\ Vu \ 2 ’ 
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s-iN P,tt _ 21 N P,tt s~iNP,tt _ 21 Nptt 

SRR2,ri d t]*\r] u \ 2 7g SRRl,i) d r)Z,\riu\ 2 ' SRR2,(VdVu) 2 76 SRR1 >( r ld , nZ) 2i 


(A.26) 


s~iNP,cc _ Q s-iNP,cc _ Q y^NP,cc _ Q s-iNP,cc _ q s~iNP,cc 

L 'SRR2,7 ldV *\r lu \ 2 ~ U ’ U SRR2,(r) dV *) 2 ~ U ’ L 'SRR2, VdV * ~ U > U SRR2,\ Vu \ 4 ~ U ’ L 'SRR2,\ Vu \ 2 


0. (A.27) 


Here x b = x c = and the functions /, (x c , Xs, x b , x t ) are defined, respectively, as 

TYIttt 772 w 


/i (x c ,x s ,x b ,x t ) 
f2 (x c ,X S ,X b ,X t ) 

h ( X c ,x s ,x b ,x t ) 
h {x c ,x s ,x b ,x t ) 
h (■ x c ,x s ,x b ,x t ) 

/e (x c ,x s ,x b ,x t ) 
/ 7 (x c ,x s ,x b ,x t ) 
h (x c ,x s ,x b ,x t ) 
h ( x c ,x s ,x b ,x t ) 


-x c x s x t 


9(x s - l) 3 (xs - 3x c ) (x s - Xt) 3 ’ 

= — 6x| + 12x| (x t + 3) — 3x| (xy (x t — 6) + 2 (x 2 + 12x t + 9)) , 

= x% (x b (x'l - 8x t + 2) + 12 (3x 2 + 9x t + 2)) , 

= x s x t ( x b (6 x t - 41) - 6 (9x t + 8)) + x t (x b (13x t + 12) + 24x t ), 
2 

X c X-f- 


9 (x t - l) 3 (x t - 3x c ) (x s ~ x t ) 3 ’ 


x s ( x b ( x t + 7x t + 2) - 6 (x t - 4) (x t - l) 2 ) , 
x s x t (12 (x t - 4) (x t - l) 2 - x b (3x 2 + 12x t + 5)) , 

X t (x b (21x 2 - 23x t + 12) - 6 (x t - 4) (x t - l) 2 x t ) , 

_ 2 x b x c x t _ 

90s - l) 2 (x t - 1 ) 2 (x s - x t ) 2 (2x c x s x t + 2x c x s + 2x c x t - x s x t ) ’ 


/io {x c , x s , x b , x t ) = x 4 s (Ax 2 - 15a:t + 6) + x| (-9x* + 17x 2 + 19x t - 12) , 

fn (x c , x s , Xb, Xt) = X 2 S (4x 4 + 2Axf - 70x 2 + 21x4 + 6) , 

fi 2 (x c , Xs, X b , Xt) = x s x t (—17x4 + 21a; 2 + 16x t - 15) + 3(x t - l) 2 x t 2 , 


/l3 (x c ,X S ,X fe ,Xt) 


2X5X 2 

9(x S - 1) 3 (X5 - Xt) 4 ’ 


fu {x c ,X S ,X b ,X 4 ) 
/l5 (Xc,X S ,X b ,Xt) 

/l6 (x c ,x s ,x b ,Xt) 
fn {x c ,x s ,x b ,xt) 


= 3x| - 6x| (x t + 3) + 3x| (x fe (x t - 3) + x 2 + 12x t + 9) , 

= — X 2 S (xb ( X 4 + IOX 4 + l) + 6 (3x 2 + 9 x 4 + 2 )) , 

= 3 x 5 X 4 (x fe (x t + 13) + 9 X 4 + 8 ) - 12x t (x b (x t + 1) + x t ), 
2 x 2 

~ 9{xt-l)\x s -x t f 


(A.28) 

(A.29) 
(A.30) 
(A.31) 

(A.32) 

(A.33) 
(A.34) 
(A.35) 
(A.36) 

(A.37) 
(A.38) 
(A.39) 

(A.40) 

(A.41) 
(A.42) 
(A.43) 

(A.44) 
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(A.45) 


/is (x c , x s , x b , x t ) = x% (x b (9 x t + 1) + 3 (x 2 - 2x t + 4) (x t - l) 2 ) , 

fi9 (x c , x s , X b , Xt) = -3 x 2 s x t (x b (9 Xt + 1) + (2a: 2 - x t + 8) (x t - l) 2 ) , 

f-20 (■ X c , x s , X b , x t ) = 3 x s x t (x b (x A t - 7x 3 + 24a: 2 - 12x t + 4) + x t (x 2 + x t - 2) 2 ) , 


/21 (x c , x s , x b , x t ) = -x A (x b (x 2 - 3x t + 12) + 9(x t - l) 2 ) , 

f -22 [x c , Xs, x b , Xt) = /is (x c , Xs, x b , Xt) + fig (x c , X S , X b , X t ) , 

/23 (x c , Xs, x b , Xt) = f-20 (■ X c , Xs, X b , X t ) + / 2i (x c , X S , X b , X t ) , 

Xt 


/24 ( X c ,X S ,X b ,Xt) = 


27 (x s - 1 f(x t - l) 3 (x s - X t ) 


3 ’ 


(A.46) 
(A.47) 
(A.48) 
(A.49) 
(A.50) 
(A.51) 


/ 2 5 (x c , xg, x b , x t ) = x b Xg (5x 3 + 2x 2 + 65x t — 12) — 2x b x 3 s (lla/ — llx 3 + 85a: 2 + 47x t — 12) , 


/26 (x c , xs, x b , x t ) = x b x 2 s (5x 3 + 130x a — 246x 3 + 574a: 2 — 91x t — 12) , 
/27 O c , Xs, X b , x t ) = -2x b x s x t (17x a - llx 3 + 49a: 2 + 95x t - 30) , 


/28 (x c , x s , x b , Xt) = +x b x 2 (~31x 3 + nox 2 - 43a: t + 24) , 

/29 (x c , x s , x b , x t ) = -180s - l) 2 Os - x t ) 2 (x t - if (x t - 4) x t , 

/30 (x c , x s , x b , x t ) = / 25 (x c , x s , x b , x t ) + / 26 (x c , x s , x b , x t ) + / 27 (x c , x s , x b , x t ), 


/31 (x c , Xs, X b , Xt) = / 28 Oc Xs, x b , Xt) + / 29 Oc Xs, x b , Xt) , 

5 x b x c xsx t (. Xg (x 2 + Xt + 2) - Xsx t (3x t + 5) + Ax 2 ) 


hi (x c ,X S ,X b ,X t ) = 
fs3 ( X c ,Xg,X b ,X t ) = 


90s - 1) Os - 3x c ) Os - x t y 

nr* ry* ^ 
tjdy (yds 


9(x t - 1) (x t - 3x c ) (x s - x t ) 


3 5 


/34 Oc £s> x b , x t ) = x 2 s (x 2 - 2o: t + 2) + a: s o: t (-3a: 2 + 6a:* - 5) + O' (3a: 2 - 6x t + 4) , 

5x b x c x s x 2 

T35 \X C , Xs, X b , Xt) = -o-o-o- 5 

180s - 1 ) 0* - 1) Os - x t ) (2x c Os 0* + 1 ) + x t ) - a: s aO 

/36 Oc ^s, X b , Xt) = x 3 g (x t - 3) + X 2 S ( -3x 2 + 5a; t + 4) , 

/37 Oc Xs, X b , x t ) = Xs {Ax 3 - I3x t + 3) + Xt (-8a: 2 + 15x t - 5) , 


(A.52) 
(A.53) 
(A.54) 
(A.55) 
(A.56) 
(A.57) 
(A.58) 

(A.59) 
(A.60) 

(A.61) 
(A.62) 

(A.63) 
(A.64) 


44 



(A.65) 


/38 ( X c ,X S ,X b ,X t ) 


Wx b x s x 2 (x 2 s (x 2 +x t + 1) - 3x s x t (x t + 1) + 3x 2 ) 
9(x s - l) 3 (x s - x t ) 4 


/ 39 ( x c ,x s ,x b ,x t ) 

/40 (Z c , 2^5) X b , *^t) 
hi (x c ,x s ,x b ,x t ) 

/42 (■ X c ,X S ,X b ,X t ) 
f43 (x c ,X S ,X b ,X t ) 


10x b x 2 (3xgx 2 — 3x 2 s x t + Xg — (x 2 — 3x t + 3) xf) 

= 9(x t — l) 4 (o;s — x t ) 4 ’ 

5x b x 2 

27 (x s - 1 ) 2 (x t - l) 3 (x s - x t ) 3 ’ 

= Xg (2x 2 — 7x t + ll) — x’g (7x\ — 19x 2 + 17x t + 19) , 

= x 2 s (llx 4 - 17x 3 t - 15x 2 t + 55 x t + 2) , 

= —xgx t (l9x 3 — 55x 2 + 53x t + 7) + x 2 (2x 2 — 7x t + ll) , 


(A.66) 
(A.67) 

(A.68) 
(A.69) 
(A.70) 
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